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I . INTRODUCTION 

The subcrltical assembly initially was utilised as an 
experimental unit to determine if a chain-reacting assembly 
were possible. As reactor technology progressed the assembly 
came to be used more and more as a test unit for proposed 
reactor designs and materials, and the majority of these test 
investigations have been conducted in the steady or equili- 
brium state. Operation of the assembly by disturbance of the 
steady state has led to new areas of investigation. 

Frequent utilization of the subcrltical assembly as a test 
unit has come about naturally. Co3t of a subcrltical assembly 
is considerably lower than that for a critical assembly since 
special safety equipment is not needed (11, p. 103). Less 
fuel and moderator are used, and experiments can be conducted 
easier in the subcrltical than in the critical assembly. 

Disturbance of the steady state may be accomplished by 
two methods, reactivity forcing and source forcing. Re- 
activity forcing can be represented as a change in the neutron 
density brought about by a small harmonic variation in the 
reactivity, whereas source forcing is produced by a variation 
in the rate at which neutrons are being added to the assembly. 

Thus far the Investigation of the dynamic response of a 
system has been limited to servomechanism theory, which 
treats the system as a "black box" . Phase and amplitude 
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response of the system In relation to a forcing function must 
be determined In order to apply servomechanism theory. The 
above relationship called the transfer function is derived by 
noting the responses as Indicated brought about by either 
reactivity or source forcing. 

It is believed however that neutron diffusion in a sub- 
critical assembly will not produce Identical phase and 
amplitude responses throughout. These variations should depend 
upon the method of varying the rate of neutrons from the 
source and the position of neutron -flux measurement. 

The purpose of this investigation was threefold. First, 
an experimental verification of the theoretical transfer 
function for a subcritical assembly was attempted using a 
varying rate of neutrons from the source. Secondly, neutron 
diffusion in the subcritical assembly that 3hows variance 
from the predicted identical phase and amplitude response 
throughout the assembly was investigated. The last area of 
Investigation was an attempt to determine assembly character- 
istics by using a sinusoidally varying source. 
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II. REVIEW 0? 'IKE LITERATURE 

A reasonable amount of literature is available on the 
development and usage of the transfer function of a critical 
assembly. Original work in the field was carried out by 
Owens, Crever and Pigott (12) in 19^3 , in which an automatic 
control unit was used to compensate for reactivity variations 
in order to maintain the desired power level in a reactor. 

Their work showed that the assembly transfer function ex- 
hibited proper characteristics to describe the reactivity 
variations in the assembly. 

Later that same year, Franz (A) developed the general 
form of the transfer function of a nuclear reactor by apply- 
ing the standard techniques of servo theory, lie applied the 
transfer function to an electronic assembly simulator to pre- 
dict assembly response to a step-function change of 
reactivity. 

The practical verification of using a subcritical assembly 
in testing reactor components was shown by Axtmann, Dessauer, 
and Parkinson (1) in 1955. A chain-reacting ’’test assembly" 
was utilized to test materials by noting effects on the kinet- 
ic behavior of the assembly. Identical tests were conducted 
in a subcritical assembly, and the results were nearly 
Identical with the subcritical assembly producing results 
faster, cheaper and safer. 
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Glass tone ( 5 ) briefly discussed the transfer function in 
a subcritical assembly . In hia development he assumed a 
sinusoidal variation in source strength and used electrical 
analogs of reactor quantities to indicate the relationship of 
electrical parameters t with the reactor parameters, reactivity 
forcing and source forcing. 

Weinberg and Schwoinler (15) developed the equations for 
the response of a critical chain-reacting assembly to a 
thermal neutron absorber which is oscillated back and forth 
inside the assembly. They further stated that at frequencies 
which are low compared to the reciprocal periods of the de- 
layed neutrons, the neutron density in the assembly rises 
and falls as a whole. However, if the frequency is much 
greater than the reciprocal period of the delayed neutrons, 
a spherical neutron wave emanates from the vicinity of the 
oscillator and is attenuated in an exponential manner with 
distance as it i© propagated. 

Weinberg and Wigner (16) extended this theory to aub- 
crltical assemblies. It was stated that the neutron-wave 
wave length is related to the material buckling. They sug- 
gested that the exponential -experiment method of determining 
material buckling could be performed by measuring the proper- 
ties of neutron waves established by a localized oscillating 
neutron source. A particularly desirable feature 
of this method is that an experiment to determine material 
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bueklins can be performed in an assembly of smaller dimensions 
than that required by the normal exponential experiment. 
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III. LIST OP SYMBOLS 



A Activity 

a Width of assembly 
a.l. Attenuation length 
B s Buckling 

(3 Fraction of average group of delayed 
neutrons 

@ Fraction of 1 th group of delayed neutrons 

C Concentration of average group of delayed 

neutron precursors 

C 0 Ena -correct len factor 

Cj, Harmonic -correct ion factor 

Cji Concentration of i til group of delayed 
neutron precursors 

D Thermal neutron diffusion coefficient 

f Thermal utilisation 

Y Inverse relaxation length 

i V^T 

K Constant 

k Effective multiplication factor 

Multiplication factor for an infinite 
assembly 

L Neutron diffusion length in assembly 

Jl Average neutron lifetime in a finite 

assembly 

J? 0 Hmm neutron lifetime in an infinite 
assembly 



counto/mln. 

In. 

in., cm. 
cm.-*® 



cm. 

in.- 1 



cm. 
sec . 
sec. 
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sec .~ 1 

sec .** x 

in. 2 , cm. 2 

neutrons/cm. 2 

radians 
sec. 

neutrons 
cm . 2 sec . 

in.*" 1 , cm.** 1 

sec . 

cm. 2 

cm ./sec . 

in. 

Subscript ; 
a Absorption 

dn Delayed neutron 
e External 

g Graphite 

1 1 th part 

mn Harmonic mn 

p Prompt 

0 Steady state 

1 Complex amplitude 
11,13,31,33 Harmonics of Y 



X Decay constant of average group of de- 
layed neutron precursors 

4- K 

X± Decay constant of the 1 * group of de- 
layed neutron precursors 

M 2 Migration area 

n Number of neutrons of thermal energies 

CJ Frequency of oscillation 

p Resonance escape probability 

/ Thermal neutron flux 

£ Macroscopic cross section 
t Time 

T Fermi age 

v Thermal neutron mean velocity 
2 Vertical distance above base of assembly 



IV. THEORETICAL ANALYSIS 



A. Transfer Function 

j 

The transfer function is defined as the ratio of the 
phase and amplitude response to a sinusoidally -varying forcing 
function. This response for a subcritic al assembly can be 
obtained from the time -dependent general diffusion equation. 
Thermal neutrons only will be considered since the experi- 
mental equipment is designed to vary the rate at which thermal 
neutrons enter the assembly. The ba3ic equation i3 (6, 

p. 101) 

,, c) n 1 c) 

v~j*- = * ~v St Kq * 1 

where $ * n v, and v is assumed to be constant. 

The source term can be separated into three parts. The 
first tern is due to the prompt neutrons resulting from 
fission. Fermi -age theory yields this part of the source 
term to be (5> p. 226) 

Sp - (1 - 0) r a 0 Q ~ Eq. 2 

where (3 is the fraction of all the delayed neutrons. 

The second part of the source term is due to the delayed 
neutrons resulting from decay of various fission products . 

A total of six groups of delayed fission neutrons are 
utilised in expressing this part of the source term. Each 
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group has Its own value of ^ ^ and C^, The net rate of 

formation of the delayed -neuetron precursors Is then 



i£i. 

at 



(3< 00 



— B 



z 



~ A I C 1 



Eq. 3 



The delayed-neutron source term using Feral -age theory can be 
expressed as (5# p. 227) 

-a- y ' 'nt 

s dn " P ® ‘ ^/_ A c l Eq, 4 

1*1 

where ZT A i is the decay rate of the precursor which is 
i**l 

numerically equal to the rate of production of the delayed 
neutrons. The Feral -age for the delayed neutrons has been 
taken equal to that of the prompt neutrons (5, p. 227). 

The last part of the source term is due to the external 
source which in this analysis supplies the source forcing 
needed. It is written as 

S e « 3 0 * S t c v M . Eq. 5 



With substitution of the above expressions of the source 
term Into Equation 1 the following equation is obtained 

D V 0 - Z ... 0 + (1- 0) k„ Z a 0 e~ tJT 

«B s t 6 . _ _ j u) t 

+ p e j- X C t + S 0 + S, e* 

i«wl 

*• -X- -ii. 

v <} t 



Eq. 6 
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The space and tine variables are considered separable 
when k remains constant (5* p. 227 ). This yields an ordinary 
differential equation for Equation 6. If the subcritical 
system is large, then v 2 0 « -B s p! (6, p. 361). 



Let 




JL 

v 



Eq. 7 



and L 2 «= D/ £ Eq. 8 

and divide Equation 6 through by £ a , to obtain the result 



- (L a B 2 + 1) 0 + (1 - (3) k^ e“ B " T pf 



+ 



p «b 2 t 
- — e 



6. 

Z 



I«1 



A i c i 




+ 



S A e * 3 U 



a 



0 ±± 

A ° d t 



Eq. 9 



With division of Equation 9 through by (1 + L 2 B 2 ) and the use 
of the relationships 



k 






e -B 2 T 



1 + L a B 2 



Eq. 10 



and 



l 



1 + L^ 2 



Eq. 11 



Equation 9 becomes 
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-pf+(l-(3)k;^*f ~T" 21 ^i 

Za 00 1«1 

3 0 4- S A ^ fc ^ d ft 
* Z a U + L*B») " d t * 

With division by i , end utilisation of 0 * n 
12 is reduced to 

, v , n v p k 

Ck(l - 0 ) - 13 - 7 *- + ~T— *Jk~ Z Ai Ci 

^ Z a K 00 l«al 

S Q + S. e" ^ a „ 

+ r a ^u + L e B e ) ” v a t • 

With rearrangement Equations 10 and 11 becosie 



-B s r 



^ k. 



~B a r T „ 

e 2 a v 



and 



i 0 - i(l + L*B S ) 



y v 
•<£ « 



or 



1 /v * Z» ^ u ♦ l* a B a ) 



Equation 13 with use of the above two substitutions 
division through by v yields 



[k(l • (3 ) - Ij “*y~ + p © ^ 
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2 T Ai < j 



i»i 



♦ S 0 + S A o 



Jwt 

ass 

d t 



Eq. 12 
5 Equation 



Eq. 13 



Eq. 14 



Eq. 15 

ana 



Eq. 16 
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In a similar manner. Equation 3 with the use of Equations 
10 and 13 is found to be 



d Ci k n 

77~ ° °1 • , E <1- w 

a t pie 



Solving for A i c i and substituting into Equation 16 yields 
[led -(3) - n -f- + P ■ 1 d Cl 



1=1 



(3i k n 



p e 






d t 



+ S Q + Sj, e 



J CJ t A n 



d t 



Eq. 18 



or 



£§- . (M , » . P £ _±|L 

d t j? 1=1 d t 



+ s 0 + Si e 



J CJ t 



Eq. 19 



since 



6 (3 k n (3 k n 

ill 7 * ' — 7 ' " 



Solving this for n gives 



n 



i 



(k-1) 



P e 



-B*r 



1 d Ci 



i«l 



d t 



d n 

TT 



H (s p + S x e* u ) 

(1-k) 



Eq. 20 
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In this particular case the value of k will remain 
constant. Neither assembly geometry nor fuel-loading changes 
will take place in the experiment . However both the neutron 



density and precursor concentration will vary as 



n « n 0 + n x e 



jut 



Eq. 21 



and 




Eq. 22 



where n x and C x are the complex amplitudes. With the aid of 
Equations IT, 21 and 22 , Equation 20 yields 



The complete solution to Equation 23 Involves using the 
six groups of delayed neutrons each with its discrete value 
of X i* (3 it and C*. The solution using an average group of 
delayed neutrons will be shown. 

Equation 23 can be separated into steady and sinusoidal 
components, and each constitutes a valid equality. The 



n 0 + n x e 





^ (S 0 

(1-k) 



Eq. 23 
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transfer function uses only the sinusoidal component, which 
with croup notation eliminated, is 



R x « 



Jwt 



l 

(k-1) 



0k 

i 



n, e 



UJ b 



- p e** sar Ac 1 e^ ^ (nj e^ ^ 

i Sj, u t 

+ . 

(l-k) 



Eq. 24 



Again using Equation IT and substituting Equations 21 
and 22 gives 

*-Lrr— li^+^eJu 4 ). X (C 0 *) 

P o i 

+ ~lk~ tC o + c i ® J U ') • Sq. 25 

The sinusoidal portion of which becoiaes 



k 0 



-E s 7* 



n, e 



UJ 



X .<■ x ® 



i cj t 



Pie 

+ 4 6J C x ^ , 



■which reduces to 



(3 k n, -a* r 

p £ e 

( A + j w) Ci 



Eq. 26 



Eq. 2? 



Substituting Equation 27 into Equation 24 yields 
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n, e 



iuJ t 



(3 k n l q v 00 A C- e*> w 



(k-1) 



(k-D 



(3 k n x 



(A + Jcj) Ci 

*S. 

a 

( l**k ) 



i jun x e J ^ * 
(k-1) 



Eq. 26 



By definition, the amplitude portion of the transfer 
function is the ratio of the response to the sinusoidally - 
varying driving source. The ratio can be expressed as 






S* 



1-k 

T~ 



Eq. 29 



or from Equation 28 





1-k 


Si 


7 



r» + — 

u (1-k) 



1 

_ A k (? 

( A tjw jll-le) 



•f 



S J cj _ 
(1-Jc) J 



Eq. 30 



Rearranging Equation 30 by putting the right side denominator 
under the common factor ( A * j£o) ( 1-k ) , multiplying through 
by x + <) k> , then dividing by A and simplifying gives the 
desired form 
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n x r 1-k 

sTl” 



2 

cj i_ + 
A(l-fc) 



1 + 



Jcj 

x 



„ CJ 

A 



El + 



Al + (3 k 

(1-k) 



Eq. 31 



B. Attenuation Length of Neutron Wave 

When the rate at which neutrons entering a subcritical 
assembly is oscillated In a sinusoidal manner at a frequency 
considerably greater than the reciprocal period of the de- 
layed neutrons a neutron wave is propagated. With an average 
delayed neutron period of approximately 10 seconds, fre- 
quencies greater than one radian/second should produce these 
neutron waves. Frequencies considerably less than the 
de 1 ayed -neutron reciprocal period allow the delayed neutrons 
to remain in equilibrium with the prompt neutrons (14, 
p. 100). However, at higher frequencies the assembly be- 
haves as if there were no delayed neutrons and the delayed 
neutrons merely contribute to the background. Hence the 
neutron wave has characteristics that are dependent upon 
the “prompt" response of the assembly. These characteristics, 
wave length, velocity, and attenuation length are all 
dependent upon the frequency of oscillation. 

Weinberg and Wigner (16, p. 437 ) developed equations 
for the three neutron-wave characteristics. The equation 
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for the characteristic that is of importance here, attenua- 
tion length, la 



a * 1 • ° 

Van s «- 32 

[«*, - l)* + ( UJJ? 0 )*W* - (kp - »ij «/• 

where u) is greater than one radian/ see ond . Attenuation length 
is defined as the distance for wave amplitude to decrease by 
a factor of e. 

If a neutron wave were propagated by using a high- 
frequency sinusoidal source oscillation then attenuation 
length could be determined. Amplitudes of the sine wave 
in selected vertical positions could be measured, and a value 
for attenuation length would then be determined from the 
equation A » K e~ s / a * 1 *, where A is the amplitude at z, and 
K is a constant. A semi-log plot of amplitude versus verti- 
cal position will yield a slope of 1/a.l. This value of 
1/&.1. is the "inverse relaxation length” of the complex 
buckling for the frequency investigated. If the phase angle 
is small at a given frequency, the attenuation length at that 
frequency would be approximately equal to relaxation length 
used to calculate material buckling. Therefore, the complex 
buckling determined by this value of attenuation length 
would be approximately equal to the material buckling. 
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V. EXPERIMENTAL EQUIPMENT 
A . Subcrltical Assembly 

The subcrltical assembly used for the experimental in- 
vestigation is shown in Figure 1 with th© east cover removed . 
Cylindrical rods of graphite measuring 7 inches in diameter 
were machined to a square from 6 inches across, leaving 
rounded comers. These blocks of AGH -grade graphite were 
then stacked in ten columns nine rows high to form the lower 
section of the assembly. For the top five rows cylindrical 
graphite rods 6-3/8 inches in diameter were machined to a 5 
inch by 6 inch cross section with rounded comers . The 
space available between the comers of adjoining blocks 
was utilized for insertion of fuel elements or measuring 
equipment . 

The graphite blocks wore assembled upon a plywood sheet 
supported by a wooden pedestal. As shown in Figure 1, the 
pedestal was divided into three accessible spaces approxi- 
mately one foot high. The two outermost spaces contained 
open-topped aluminum tanks, filled with water, which ex- 
tended the width of the assembly. They served to moderate 
the neutrons and to reduce radiations coming directly from 
the sources located in the oscillator unit beneath the 
assembly. Figure 2 shows the oscillator unit and selective- 
count device which was placed in the center access section 
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Figure l. The suberitical assembly 
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Figure 2. The oscillator unit and selective- 
count device 
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beneath the assembly. 

All four sides and top of the assembled graphite were 
covered with 10-mil cad.nl am sandwiched between a G. 375-inch 
thickness of plywood and 0.125 Inches of masonite. Approxi- 
mately 95 percent of all thermal neutrons were stopped by 
this cadmium sheet, giving a “black boundary" to thermal 
neutrons. A wooden framework braced the north and south 
sides of the assembly and both the east and west faces wore 
removable for easy access. 

Natural uranium in the form of cylindrical rode 1.0 
Inches In diameter and o.Q inches long served as the fuel. 
These rods were encased by 2S aluminum cans with a O.OtO- 
inch wall thickness and end caps 0.200 inches thick. Each 
can of uranium was then helically wound with a spacing wire 
and placed into a 6lS aluminum tube, 62 inches long with an 
outside diameter of 1.375 inches and a wall thickness of 
0.035 inches. Approximately 10 feet of the 2S -aluminum 
spacing wire was needed for each tube, A total of seven 
canned uranium rods were placed in each tube to constitute 
one fuel element. With the fuel elements placed in every 
other hole as shown in Figure 1 an 8. 5-inch square lattice 
in the lower region of the assembly was formed. This lattice 
configuration was maintained throughout the experimental 
investigation. 
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B. Oscillator Unit 

The oscillator unit and selective-count device is shown 
in Figure 2. The drum-typo assembly at the right end is 
the oscillator unit itself and is shown schematically in 
Figure 3. This unit was designed to contain the five one- 
curie Pu-Be sources and to vary the rate at which thermal 
neutrons were emitted from the source in a sinusoidal 
manner. 

Each cylindrical source measuring 1 inch in diameter 
and 1-3/B inches high was placed in the source retainer. 

A non-rc tat ing -paraffin Jacket surrounding the retainer was 
constructed of 1/16-inch aluminum sheet rolled into a 
diameter of 7-3/3 inches. End capo of 1/0-inch aluminum 
were attached, giving an overall length of 14.5 inches. 
Twenty-five pounds of white paraffin wore poured into the 
Jacket through a filler cap, and precautions were taken to 
reduce the possibility of shrinkage voids, Moderation to 
supply the thermal neutrons do aired was provided by the 
three inches of paraffin. The paraffin Jacket and the bear- 
ings for the pattern cylinder were both supported by the 
source retainer which extended to the support brackets. 

An outer rotating cylinder, the pattern cylinder, wa3 
also formed from 1/8-inch aluminum sheet. Its dimensions 
were an outside diameter of 8 Inches and a length of 15 
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ROTATING CYLINDER 
WITH PATTERN 



SOURCE RETAINER 



NONROTATING 
PARAFFIN 
JACKET 



PARAFFIN 




SUPPORT 

BRACKET 



( 360 °) 



Figure 3. Oscillator unit and cadmium -pat tern detail 



inches. End plates made from 1/4-inch aluminum were fitted 
and screwed to a collar fitted on its bearing outer-run 
casing. The inner-run casing of the tearing rested on the 
source retainer for support. To give easy access to the 
inner section of the pattern cylinder the end plates were 
attached by screws. 

A sine pattern with the dimensions as shown in Figure 3 
was cut from 10-mil cadmium. Three thicknesses, giving a 
total of 0.030 inches, were cemented to the inner surface of 
the pattern cylinder. Most of the thermal neutrons were 
stopped and a small fraction of the fast attenuated by this 
thickness of cadmium (10, p. 12), 

Point A of the pattern was the wide-open position, allow- 
ing the maximum number of thermal neutrons to escape from the 
oscillator unit. The closed position, point B, stopped moat 
of the thermal neutrons. When point A was at top-dead 
center, the oscillator unit was at the aero-degree position. 
The loO -degree position occurred when point B was at top-dead 
center. lienee the exact position of the sine pattern of 
the oscillator unit was designated by a given angle from 0 
to 360 oscillator degrees. 

Measurements taken within 1.5 inches of the oscillator 
unit at the two extreme positions described above showed 
that the paraffin Jacket moderated a satisfactory percentage 
of neutrons to be used in the experimental investigation. 
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The value obtained at position A minus that at B was 58 per- 
cent of the value at position A. 

C . Counting Apparatus 

In order to measure the neutron flux during a designated 
number of oscillator-unit degrees a selective-count device 
was constructed. A positive, non-slip linkage had to be 
supplied from the oscillator unit to a position removed from 
the neutron sources. The equipment needed to select the 
number of oscillator degrees in which counting would take 
place was located at this "removed" position. 

A bicycle sprocket was bolted to the end plate of the 
pattern cylinder. A 110-inch bicycle chain was used to 
connect this sprocket on the oscillator unit to a sprocket 
of the same diameter mounted on a countershaft at the removed 
position. Position by degrees of the oscillator unit was 
then transmitted unchanged to the "removed" position where 
work could be performed without the intense radiation. 

The countershaft was a polishing -head spindle which was 
mounted on a 3/4-inch plywood extension connecting the 
oscillator unit to the remainder of the equipment. At one 
end of the spindle shaft a 6-inch V-belt pulley was locked 
to the shaft, and the bicycle sprocket was in turn bolted 
to this pulley. At the opposite end of the shaft was locked 
a plywood plate, 6-3/4 inches in diameter. The polishing- 
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head spindle was driven by a 1/4-inch Black and Becker 
electric drill connected to it by a leather belt. Rotation 
frequency was controlled by varying the voltage to the drill. 
In order to improve the frequency control at alow speeds 
(below 50 RFM) a second polishing-head spindle was interposed 
as shown in Figure 2. These two configurations then gave a 
speed range for the oscillator unit of from 4 to 1000 RFM. 

A 1/4-inch plywood plate of 8. 14-inches diameter with a 
30-degree cam hump wt fabricated. To prevent wear a metal 
atrip was cemented to the cam. A similar 10-degree cam was 
also constructed* The cam, with a 1/2-inch centering hole, 
was held securely by a nut against the plywood plate on the 
spindle shaft. A compass card oriented to the exact 
oscillator degree was secured to this plywood plate, the 
cam was set at any desired oscillator position by loosening 
the cam-holding nut. 

As shown in Figure 2 a microswitch was positioned adja- 
cent to the cam. The switch was wired in the normally -open 
position to the preset-time relay of a Radiation Instrument 
Development Laboratories, model 206, scaler. The cam was 
positioned by use of the compass card to make contact with 
the switch at a specified oscillator degree. When the sine 
pattern of the oscillator unit reached this degree the cam 
made contact and the switch was closed. This closed the 
preset-time relay allowing the scaler to begin to count. 
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When the microawitch broke contact the relay of the scaler 
opened, ceasing the count* 

A Radiation Counter Laboratories boron trifloride pro- 
portional counter, model 10 ^ 03 , was inserted midway into the 
assembly through holes cut into the east cover. It in turn 
was connected to the Radiation Instrument Development 
laboratories scaler. 

The oscillator unit and selective -count device was in- 
serted from the south end of the assembly until the oscillator 
unit was at the center. As shown in Figure 4 this placed 
only the oscillator unit itself beneath the assembly; the 
remaining equipment was outside to facilitate changing the 
can setting. 
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Figure 4. Experimental set up 
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VI, EXPERIMENTAL PROCEDURE 
A. Description of Typical Run 

She fuel elements were loaded into the assembly as shown 
in Figure 1, Measure “onto were taken only In the center of 
the assembly. Holes drilled on the east-cover centerline 
matched assembly holes between the graphite blocks. The 
holes were one foot, two feet, three feet, four feet and 
four feet eleven inches above the base. Hereafter these 
holes will be referred to as holes one through five 
respectively. 

Tho neutron detector was inserted midway into the assembly 
and the oscillator unit set at 0 degrees. With the 30- 
degree cam set at a mm position of 0 degrees, initial 
contact with the mlcroswltcb was made at 3^5 degrees. 

Oscillator speed was adjusted to the desired RPM by re- 
ferring to a stop watch and counting revolutions. When de- 
sired speed was obtained the watch was started and the 
preset-time switch on the scaler was actuated. Runs varied 
in duration from one to five minutes depending on position 
and oscillator speed. A run was completed when a designated 
number of counting periods had elapsed rather than ending at 
an exact time interval. For example, a run for 3 minutes 
at 18 BPM requires 5*4 counting periods; when 84 periods had 
elapsed the preset-count switch was turned off even though 
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the time my have beer* a few seconds cn oither side of 3 
minutes. It was difficult to control the speed exactly and 
this procedure eliminated varying numbers of counting periods 
in a circle scries of runs. 

The counts were recorded from the scaler} the drill 
turned off; the cam reset to a mean position of 30 degrees 
and the entire counting procedure repeated. The counting 
procedures continued in 30 -degree increments until the entire 
360 degrees had been surveyed yielding one single series of 
runs* 

Speeds (frequency of oscillation) used in the investiga- 
tion were selected at regular intervals on a logarithmic 
scale. Huns were attempted throughout the entire range of 
the equipment . The lowest constant speed attainable was 
l t RPM. An upper limit on speed was Imposed by the cam and 
microswitch arrangement, for in speeds beyond ISO RBI the 
switch hammered the cam es the switch spring forced the idler 
back: to the surface of the cam. As runs were made in the 
higher hole positions at high speeds the data became unusable. 
This was caused by the exponential attenuation of neutron 
flux up the assembly end a consequently smaller difference 
fro® the highest to the lowest count rate in a series of 
runs. Vfhen this occurred the counting rates when plotted 
against degrees did not yield a satisfactory distribution. 

Original data for the various runs completed is shown in 

Tables 3 through 12 in the Appendix. 
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B. Determination of Phase Angle 

The count rates for a single sex'les of runs were plotted 
on graph paper. Standard deviation for each point was also 
plotted to facilitate drawing a representative curve through 
the points. All the plots that showed a trend in distribu- 
tion gave a sinusoidal result. Figure 5 shows a plot of the 
data for all holes at 2. 51 radlana/cec ond , and as mentioned 
previously some unusable data was taken and the basis for 
declaring it unusable was large scatter in the plotted 
curves . 

From the general trend of the plot and especially from 
the position of the highest and lowest counting rates the 
phase angle for that series of runs was determined. To ensure 
that an accurate phase angle was being obtained a 10-degree 
can was utilised. Five-minute counts were taken every ten 
degrees in the region of the maximum counting rate, however 
use of the 10-degree cam was restricted to the lower speeds. 

It was assumed in determining the phase angle that the 
oscillator unit when in the O-degree position during rota- 
tion and the counting equipment would yield the highest 
counting rate. To validate this assumption an entire series 
of calibration runs with the neutron detector held 1.5 
Inches from the oscillator unit in position beneath the 
assembly was made* Table 13 in the Appendix and Figure 5 
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indicate the results which show that there was no "built in" 
phase angle. 



C. Determination of X by Amplitude Ratios 

The inverse relaxation length for the thermal -neutron 
flux in the assembly, X , was determined by amplitude ratios 
for various frequencies of oscillation. The value actually 
obtained was the inverse attenuation length, but as previously 
noted, under certain conditions these values are approximately 
equal. The experimental conditions were such that the above 
statement is true, and will be demonstrated in the Discussion 
of Results. 

Uncorrected -amplitude ratios were used in place of foil 

- 2T 2 

activities in the equation A « Ce . With the method of 
least squares an Initial uncorrected value for X was found 
(2, p. 166). Both end-correction and harmonic -correction 
factors* were determined and applied to the amplitude ratios 
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with this initial value of x . The origin of the Z-axia was 
taken at tiie base of the assembly; therefore hole one was at 
a position of z « 12 inches. A second value for x was 
determined with the corrected amplitude ratios . This process 
was repeated until the value of X used in the correction 
factors agreed with the X found by using the corrected 
amplitude ratios. Computations showed that harmonic cor- 
rections were significant only for the three lower holes and 
end corrections for the two upper. 

It was believed that harmonic corrections were unnecessary; 
therefore a value for X was determined in a similar manner 
using the end corrections only. 

D. Determination of X by Median -flux Levels 

Inverse relaxation length was also determined by using 
the median-flux level. Uncorrected median -flux levels wore 
used in place of foil activities in the equation A * Ce" 2 . 

A procedure was followed as explained in the preceding 
section to obtain a value for X . Values were determined 
by using only the end-correction factors. 

E. Determination of Amplitude Ratio 

Basically, the need of the amplitude ratio was to deter- 
mine if the amplitude of the neutron sine wave varied with 
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the frequency of the source oscillation. To eliminate any 
’’built-in” variation due to the equipment, the calibration 
runs next to the oscillator unit were used as the basis for 
determining the amplitude ratios. The amplitude ratio was 
defined as the amplitude of the wave obtained for a particu- 
lar position and speed divided by the calibration-run 
amplitude at the same speed. Both values of amplitude were 
obtained from the plots of the data, and to standardise the 
results, all amplitudes were placed on a one-minute basis. 

After the amplitude ratios were obtained from the data, 
an end -correct ion factor was applied. Press the iterative 
process explained previously, a value for the end -correction 
factor was obtained. No end-correction factors were applied 
to the calibration runs. The end -correction factors were 
divided into the amplitudes obtained in the hole positions, 
and it was found that factors were significant only in the 
top two holes. 

There is a definite difference between the theoretical 
amplitude ratio obtained from Equation 31 and the amplitude 
ratio obtained from the experiment. Theoretical amplitude 
ratio was defined as the ratio of amplitude of the observed 
run divided by that of an infinitely slow run, whereas the 
experimentally-determined amplitude ratio is the amplitude 
of the observed run divided by a calibration run. Therefore, 
numerical values of the two ratios will not agree, but they 
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will agree in showing any change in amplitude ratio with fre- 
quency . 



F. Determination of Median-flux Level 

The median-flux level is the average flux at a particu- 
lar hole position and oscillator frequency. From the plots 
of original data, the amplitude and the lowest value of 
count rate found in the graph were obtained. The median -flux 
level was then taken as half the amplitude added to the 
lowest count rate. End-correction factors were applied to 
the amplitudes for the top two holes. 
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VII. RESULTS 

The original data for all runs completed including the 
calibration is given in Tables 5 through 13 in the Appendix* 
Figure 5 shows a representative plot of original data at a 
frequency of 2.51 radians/second for all hole positions and 
calibration. 

Phase angle was determined for each run and the results 
are as given in Table 1. All the calibration runs proved 
that there is no "built in" phase angle in the equipment. No 
phase angle is observed in any hole position above 0,63 
radians/second. Approximately a 20-degree phase angle is 
found to exist in all five hole positions at 0.63 radians/ 
second. For the only hole position investigated, number two, 
a 30-degree phase angle is found at 0.42 radians/second. 
Verification of the phase angles in hole two up through 4.19 
radians/second was obtained by use of the 10-degree cam. No 
change of phase angle is observed in proceeding to higher 
positions while taking measurements at the same oscillator 
frequency . 

Amplitude ratios were obtained and the results are as 
given in Table 2. Figure 6 is a plot of measured amplitude 
ratio versus oscillator frequency for hole positions one, 
three and five. Amplitude ratios do not show a definite 
trend with frequency, but do remain approximately constant 
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Table 1. Phase angle (degrees) 



Frequency 

radians/second 




Hole position 






Calibra- 

tion 


One 


Two 


Two a 


Three 


Four 


Five 


0.42 


«»«• 


30 


30 








0 


0.63 


20 


20 


20 


20 


20 


20 


0 


1.03 


0 


0 


0 


0 


0 


0 


0 


1.68 


0 


C 


0 


0 


0 


0 


0 


2.51 


0 


0 


0 


0 


0 


0 


0 


4.13 


0 


0 


0 


0 


0 


0 


0 


6.28 


0 


0 


m*m 


0 


— 


m •• 


0 


7.33 


0 


0 










0 


9.42 


0 


0 










0 


12.55 


0 


0 










0 


18.85 


0 












0 


a Runs eccaupllahed with 10 -agree 


cam* 







Table 2. Amplitude ratio 


Frequency 

radians/second 






Hole position 




One 


Two 


Three 


Four 


Five 


0.63 


0.0710 


0.0330 


0.0128 


0.00610 


0.00418 


1.05 


0.0678 


0.0440 


0.0139 


O.OO063 


0.00279 


1.68 


0.0780 


0.0243 


0.0083 


0.00450 


0.00216 


2.51 


0.0810 


0.0180 


0.0090 


0.00492 


0.00274 


4.19 


0.0690 


0.0315 


0.0109 


0.00454 


0.00252 


6.28 


0.0550 


0.0259 


O.OO83 


-- 




7.33 


0.0530 


mm 


mm 


-- 


4mm+ 


9.42 


0.0673 


0.0144 


mm 


mm 


mm 


12.55 


0.0678 


0.0218 


mm 


m m 


m mm 


18.85 


0.0718 


mm 


mm 


mm 


mm 



AMPLITUDE RATIO 
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FREQUENCY- RADIANS/SECOND 



Figure 6. Amplitude ratio vs. frequency 
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at each hole position. A high coefficient of variation is 
evident for values above hole one. 

Amplitude ratios were also plotted against vertical 
distance in the assembly. The five frequencies investigated 
in all five hole positions are represented on this plot. 
Figure 7. The line slopes drawn correspond to the ar values 
obtained by the method of least squares for the set of 
points used. 

A median-flux level for each hole position and frequency 
was determined. Table 3 gives the results of this computa- 
tion, which also includes the calibration runs. To demon- 
strate the change in median -flux level with oscillator fre- 
quency a plot of these two variables was constructed for the 
first three hole positions and calibration run. Figure 8 
shows that the slopes on logarithmic paper of these lines 
for the four positions are equal. 

Table 4 gives the various values of <f obtained by using 
the amplitude ratios. Uncorrected values gave an average 
of 0.0710 Inches'’ 1 j end -and -harmonic -corrected values gave 
an average of O.063I inches” 1 , and end-corrected values gave 
an average of O.G7G2 inches” 1 . End-corrected y values were 
plotted versus frequency in Figure 9. Values for r arc 
approximately constant with frequency . Large coefficients 
of variation are noted for all the Y values plotted in 
Figure 9* 
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Table 3. Median-flux level (total counts) 



Frequency. 

radians/ 

second 




Hole position 




Calibra- 
tion 
(1 min.) 


One 

(2 min. ) 


Two 

(2 rain. ) 


Three 
(3 rain.) 


Four 
(4 min. ) 


Five 
(5 niin. ) 


0.63 


11,534 


5,410 


3,535 


2,105 


1,150 


13,585 


1.05 


11,275 


5,635 


3,730 


2,115 


1,265 


13,425 


1.63 


11,350 


5,270 


3,473 


1,99S 


1,085 


13,245 


2.51 


11,150 


5,315 


3,590 


2,098 


1,163 


12,947 


4.19 


10,857 


5,113 


3,450 


2,180 


1 , 09 S 


12,730 


6.28 


10,980 


5,100 


3,500 




•am 


12,536 


7.33 


10,634 


•* mm 




— 


tarn m* 


12,495 


9-42 


10,600 


5,078 


•m mm 


— 


«p- 


12,339 


12.55 


10,490 


4,8'30 


mmw 


— 




12,264 


18.85 


10,400 


— 


mm «» 


— 


— 


11,935 



Table 4. 


Inverse -relaxation length (inches) - 


1 




Method 

of 

computa- 

tion 


Frequency-radlans/second 






0.63 


1.05 


1.68 


2.51 


4.19 


Average 


No cor- 
rections 


0.0642 


0.0714 


0.0764 


0.0690 


0.0739 


0.0710 


End cor- 
rections 


0.0640 


0.0704 


0.0759 


0.0676 


0.0729 


0.0702 


End and 
harmonic 
cor- 
rections 


0.0582 


O.O 638 


C .0690 


0.0626 


0.0664 


O.C 631 



Figure 7. Amplitude ratio vs. vertical distance 
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Figure 8. Median-flux level vs. frequency 
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Figure 9. Inverse relaxation length vs. frequency 
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A value for ^ of 0.0681 inches" 1 was obtained by using 
the rjsthod of median -flux levels. 



* 
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VIII. DISCUSSICM OP RESULTS 

It was the original belief that a phase angle would be 
observed at any oscillator frequency above 0.6 rad iano/sec ona 
(13, p. 33). After a aeries of rune had been completed, this 
was found not to be true. Previously, the values of k and J? 
had not been known for the assembly used, and had to be 
assumed. This led to the erroneous prediction. 

The theoretical analysis was performed using one average 
group of delayed neutrons. Figures 10 and 11 show the 

theoretical phase angle and amplitude ratio respectively as 

* 

determined from Equation 31. A value for k of 0.63 was used. 
With f equal to 0.8 7 (9, P* 56), and B 2 equal to TO x 
10 cm® (9, p. 5Tb) a value of j? was determined to be 
O.OOO98 seconds from the equations 

L* « I^* J (l-f ) and ^ « 1 4.°L i *’B SE 

end an assumed value of i 0 * 0.001 seconds. The total 
fraction (3 of the delayed neutrons was set equal to 0.0064 
(T) with a corresponding value of A equal to 0.08 seconds"*. 

It can readily be seen from Figure ID that a phase angle 
would not be observed in the assembly until an oscillator 
frequency of 40 radians/second is reached. A phase angle 

*Beck, D. M. , Ames, Iowa. Data from partially completed 
M. S. Thesis. Private communication. 1959. 
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10. Theoretical phase angle vs. frequency 
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Figure 11. Theoretical amplitude ratio vs. frequency 
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does not appeal’ for this investigation until very high fre- 
quencies because the value of k is low. If k were 0.9 a 
phase angle would be observed at a frequency of 20 radians/ 
second. The overall shape of the theoretical curve remains 
the sane with varying values of k. 

Huns were completed to the lower and upper frequency 
limits of the oscillator and counting equipment . The upper 
limit (19 radians/second) is short of the required frequency 
for an observable phase angle. In complete agreement with 
theory, no phase angle was observed in any of the hole posi- 
tions above 0.63 radians/seecnd. Use of the 10-degree cam 
verified these results. 

Readings were taken in all hole positions for a frequency 
of O .63 radians/second . A consistent phase angle of approxi- 
mately 20 degrees is found. At 0.42 radians/second the phase 
angle In hole two is approximately 30 degrees. The theory 
advanced does not predict a phase angle at these low values 
of frequency. Further investigations were not conducted into 
this low frequency range since the equipment did not provide 
adequate speed control below 4 RFM. 

Theory makes no distinction of position in relation to 
phase angle; therefore the phase angle should be identical 
in all areas of the assembly. Ho change in phase angle with 
hole position was observed during the investigation. Surveys 
were taken only up the center of the assembly so the 
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proof of the constancy of phase angle is limited to this con- 
dition. With the method of source forcing used it is believed 
that some large variations would be found at the edge of the 
assembly near the base. Simply noting that the pattern 
position facing these extreme positions will vary as much 
as degrees from the recorded position indicates the 
physical reason for the nonconformance. 

As shown in Figure 5* the curve of the neutron flux 
measured at the various hole positions io shown to be 
sinusoidal. Scatter of points increases as higher hole 
positions are investigated. The range of frequency for 
usable data narrows as measurements are taken higher in the 
assembly. Statistical variations, and inaccuracies in the 
detector and scaler distort the counts received. Due to the 
already small amplitude of the sine wave, any distortions 
cause a scattering of points that make the readings unusable. 
This situation is aggravated by higher frequencies. As the 
oscillator frequency increases, the neutron waves are propa- 
gated quite close together. Any disturbance or 
irregularity in the entire assembly or counting system will 
cause an incorrect measurement to be taken. For example, a 
measurement is to be taken at a specified oscillator-unit 
degree. A small disturbance or irregularity will cause the 
neutron wave to be delayed in reaching the detector. At a 
high frequency this delay is significant since the portion 
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of the cine wave actually surveyed nay be 30 to 40 degrees 
removed from the portion desired. 

As shown in Figure 5 » the standard deviation at hole on© 
is reasonable and acceptable. However as readings are taken 
at higher positions in the assembly the coefficient of varia- 
tion caused by standard deviation Increases even though 
longer counting times are allowed. This is caused by lower 
count rates due to the exponential attenuation of the thermal- 
neutron flux. An average coefficient of variation for all 
the runs completed ranges from one to three percent in going 
from hole one to five. It is seen then that to obtain 
coefficients of variation of less than two percent at the 
higher hoi© positions, very long counting times are needed, 
in the order of 10 minutes at hole five. 

Amplitude ratios were determined and are given in Table 
2, These amplitude ratios for holes one, three and five were 
plotted against frequency in Figure 6. A least squares analysis 
was applied to the data for holes three and five, and it was 
found that the line Joining hole-three data had a slope of 
-0.184 +0.225, and the line joining hole-five data had a slope 
of -0,231 40.344. This indicates that a horizontal line 
representing a constant value of amplitude ratio with frequency 
lies easily within the standard deviation for the slopes found 
for both holes. It Is believed that the amplitude -ratio values 
are approximately constant as predicted by Figure 11 and that 
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the slopes shown on Figure 6 are caused by the high standard 
deviation rather than variance of amplitude ratio with fre- 
quency* 

Figure 11 has the overall shape that the amplitude ratio 
versus frequency would have for any value of k* For a k 
greater than 0.53, an amplitude ratio less than one would be 
noted at frequencies less than indicated on Figure ll. 

Scatter of tho points in Figure 6 is due to a very high 
coefficient of variation. As was mentioned previously, some 
scatter of points la present at all positions especially in 
the higher hole positions at higher frequencies. Error occurs 
when fairing in the curves to fit the data. A very signifi- 
cant error occurs when the standard deviation of the points 
is large compared to the amplitude. For example , at hole 
three a standard deviation of approximately 60 counts is 
foundj however an amplitude of only 300 counts exists. 
Coefficient of variation caused by standard deviation of the 
amplitude is then approximately £0 percent. The condition Is 
exaggerated at higher hole positions such that a 36 to 42 
percent coefficient of variation in hole five Is noted. 
Significant errors appear for the amplitude ratios above 
hole one. Considerably longer counting times are needed at 
all higher hole positions to give satisfactory results. 

A counting time of approximately 10 minutes would give a 13 
percent coefficient of variation for hole four. 
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Even though the coefficients of variation for the ampli- 
tude ratios are quite high, the results are believed to be 
sufficiently accurate and repetitive to support the conclu- 
sions drawn. 

Values of the median-flux level are given in Table 3 end 
the plot of level versus frequency is shown as Figure 8. 
Deviation was so slight In holes four and five that they were 
not included in the figure. 

The transfer-function development does not predict a 
decrease in the median flux, yet all four curves show that 
median-flux level does decline as frequency increases. Slopes 
on the logarithmic plot (Figure 8) of all four lines are 
found to be equal. If each line is represented by the equa- 
tion jZf ® K. CO n , where ft is a constant and n is the slope of 
the line, then the rate of neutron-flux decline with fre- 
quency is 

<X$ Kn co _ n0 

A cj " 

The relative rate of decline is 

6,0 n 
-q.^- « CO 

9 

With identical values of n, and th® same values of CO , each 
hole position has the same relative rate of decline. Th© 
decline is therefore not a function of the assembly. 
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The decline is attributed to either the method of source 
forcing or to the counting equipment. No logical reason can 
be advanced for the unit itself to cause the decline, but 
the counting equipment could easily have caused this un- 
predicted result. The number of activations of the preset- 
count relay la in proportion to the frequency. A small relay 
lag would then become important at higher frequencies and 
significantly reduce the amount of time the scaler was count- 
ing. Furthermore, investigation shows that the reduction of 
thermal neutrons and fast neutrons counted is in proportion 
to their total count©. The deviation from theory la there- 
fore said to be a mechanical irregularity of the counting 
equipment rather than a function of the assembly. 

Amplitude ratio© were not affected by the decline of 
count© with frequency since the calibration run© were used 
a© the denominator of the ratio. 

When plotted on semi-log paper, amplitude ratio versus 
vertical distance gives a reasonably straight line. Figure 7 
©hows the result© for frequencies from 0,63 to 4.19 radians/ 
second. A© noted previously, the large coefficient© of 
variation of data for the higher hole positions result in 
some scatter of the amplitude ratios. It is seen, however, 
that the points nearly all lie close to the straight line© 
indicating that the point© can be considered reliable. 
Standard deviation of the amplitude ratios is indicated on 
the figure. 
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Values for inverse relaxation length, X , were computed 
for the frequencies. These values obtained by an Iterative 
process are listed in Table and it is found that values 
range from 0.0582 to 0.C759 inches** 1 . A value for X was 
found without using either the end or harmonic corrections} 
the average value was 0.0710 inches** 1 . When both corrections 
were applied, an average value of O.O631 inches** 1 was found. 
It is believed however that the higher harmonics have been 
greatly lessened with this particular arrangement of sources 
in the oscillator. The sources are placed in a line and are 
surrounded by three inches of paraffin which thermallzes 
many of the neutrons produced. The paraffin, lined sources, 
and two water-filled tanks in the outermost pedestal spaces 
help to scatter the neutrons more uniformly beneath the base, 
so that a plane source of thermal neutrons is approximated. 
Dope hie , Leonard, Neve de Mevergnles and Tavernier (3) 
showed in their work with an exponential assembly that dis- 
tributing four sources in a square form reduced the harmonic - 
correction terra to 5 percent or one-fourth the correction 
needed when the sources were placed together. 

With application of just the end corrections an average 
value for X of 0.0702 inches** 1 was determined. Hayes (8, 
p. 27) found by the method of foil irradiation that the 
assembly in this configuration had a X of 0.0705 inches" 1 . 
Harmonic -correction factors would appear then not to be 
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needed, or if used the full correction should not be applied. 
Line slopes on Figure 7 are the ¥ values determined by using 
only the end corrections. 

Values of inverse relaxation length obtained by using 
end corrections only were plotted with frequency on Figure 9. 
A standard deviation for ¥ was obtained by using the sethod 
of propagation of precision indexes as presented by Worthing 
and Geffner (17# P* 209) and are indicated on the figure. 
Large coefficients of variation are again noted for ¥ . 

This is due to the high coefficients of variation of the 
amplitude ratios, especially at the higher hole positions. 

By using longer counting times the relative standard devia- 
tion of the amplitude ratios would be greatly reduced and 
statistically better values of ¥ would follow. For example# 
if an average relative standard deviation of 5 percent were 
found for the amplitude ratios, then the coefficient of 
variation for ¥ in this experiment would decrease by a 
factor of approximately six. 

It is seen from Figure 9 that there is no trend in ¥ 
with increasing frequency. A constant value of ¥ should be 
found, but due to the high standard deviations values differ- 
ing slightly from the average of 0.0702 inches"' 1 are noted. 

To substantiate the values of ¥ determined by using 
amplitude ratios an average value of ¥ was found by using 
the median-flux levels. As shown on figure 12 the values of 



MEDIAN-FLUX LEVEL - COUNTS/MINUTE 
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VERTICAL DISTANCE - INCHES 



Figure 12. Median-flux level vs. vertical distance 
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the median flux are grouped close together. An average value 
was obtained from each group and a value of X computed from 
them by using end corrections only. Standard deviation of 
the average points is approximately the size of the symbols, 
therefore the value of Y obtained is reliable. From Figure 
12 a value of Y was found to be 0.0681 inches -1 in compari- 
son to the average value of 0.0702 inches"* 1 computed by using 
amplitude ratios. Agreement within 3 percent between the 
two methods indicates who amplitude-ratio method to be 
acceptable . 

The validity of the statements concerning the attenua- 
tion length and 1/y was investigated. Values used in the 
computation were 'r « 350 cm 3 (16, p. 331)# ftp * k(l « (3 ) » 
0,526, L 3 * 325 cm 3 , i o * 0.001 seconds and co «* 20 
radians/second. A value was determined from Equation 32 for 
attenuation length and found to be 14.8 inches. The inverse 
of this number is 0.06/6 inches -1 which agrees favorably 
with the values of 0.0702 inches -1 and 0.0681 inches'* 1 ob- 
tained by the methods of amplitude ratios and median-flux 
levels respectively. 

Frequencies utilised in the experiment were small enough 
to cause the attenuation length to be approximately equal 
to the relaxation length. From Equation 32 it is noted that 
a frequency of 20 radians/second produces negliglbe change 
in the value of attenuation length as compared to very low 
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frequencies. Hence from the previous computations ana dis- 
cussion it appears that the attenuation length and conse- 
quently the complex buckling for the frequencies investigated 
are approximately the same as relaxation length ana material 
buckling. Therefore it would be proper in this case to 
utilize the values of attenuation length for determining S . 
Even with the combination of varying and steady neutron flux, 
computation of S by ©amplitude ratio and median-flux level 
is feasible since the frequency never is large enough to 
produce significant changes in S * 

For assemblies with much larger values of k p , frequencies 
used in this experiment would produce changes in the attenua- 
tion length and give varying values for the complex buckling, 
therefor*© the inverse attenuation length would not be equal 
to the inverse relaxation length. 
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DC . CONCLUSIONS 

A theoretical and experimental investigation into neutron 
diffusion in the subcritieal assembly with sinusoidal source 
oscillation was conducted at Iowa State College. An oscilla- 
tor unit produced the sinusoidal source oscillation of the 
thermal neutrons, while & selective -count device allowed 
discrete sections of the sine wave to be observed and counted. 
Neutron-diffusion information was extracted from plots of 
these counts. Conclusions drawn from the Information are that; 

1. No phase angle is observed for oscillator frequencies 
between 0. 63 and 18.85 radians/second. This is in agreement 
with theory. 

2. An unpredicted phase angle of 20 degrees is observed 
at a frequency of 0.42 radians/second and an angle of 30 
degrees at 0.63 radians/second. 

3. There Is no phase-angle change at higher hole posi- 
tions on the centerline of the assembly at the same oscillator 
frequency, which ranged from 0.63 radians/second to 6.28 
radians/second. This is in agreement with theory. 

4. Amplitude ratios for a given hole are approximately 
constant for varying frequencies from 0.63 radians/second 

to 6.28 radians/second. This is in agreement with theory. 

5. The median-flux level decreases as frequency in- 
creases. Mechanical operation of the counting equipment 
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causes this change rather than any assembly peculiarities or 
characteristics. A change in level is not predicted by 
theory. 

6. The oscillator unit thermal iscs enough neutrons to 
give satisfactory readings up through hole five up to 4.19 
radians/oecond. Any frequency readings above 12.55 radians/ 
second are confined to hole one. 

7. An inverse relaxation length of 0.0702 inches" 1 
was determined by using amplitude ratios. 

8. An inverse relaxation length of 0.0681 inches" 1 wau 
determined by using median -flux levels. 

9. The attenuation length of the neutron wave is equal 
to a constant value of 1/y for the assembly and frequencies 
investigated. 

10. Longer count times are needed to reduce coefficients 
of variation in the amplitude ratios. In all other respects 
the experimental equipment and procedure will yield satis- 
factory results within certain frequency limits for phase 
angle, amplitude ratio and inverse relaxation length. 
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X. SUGGESTIONS POE FURTHER STUDY 

From a review of the literature It is obvious that neutron 
diffusion in a eubcritical assembly with sinusoidal source 
oscillation has not been fully investigated. 

The equipment used in this experiment can be improved in 
several ways. Another speed reducer could be introduced to 
allow speeds below 4 RFM to be utilized. The cars could be 
replaced by a phototube device that would trigs®** the scaler 
at the higher frequencies. An automatic revolution counter 
would be necessary to measure the frequency at the higher 
speeds. A gating circuit would be more desirable to inter- 
rupt pulses going to the scaler than using the preset-time 
relay which Is burdened at high frequencies. Still another 
possibility might be to us© an ionisation chamber connected 
to a sensitive electrometer. The wave produced by a re- 
corder connected to the electrometer could then be analysed 
for data desired. 

In this experiment, only the center of the assembly was 
investigated. Measurements could be taken throughout the 
assembly to see if the thermal -neutron flux behaved in the 
same manner as th® flux produced by a steady source* Further 
theoretical and experiments’ll work with neutron waves could be 
accomplished with an investigation into the change in prompt 
buckling by source forcing. Speeds considerably below and 
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above those Investigated here could also be utilized to ex- 
pand the experimental verification of the theoretical transfer 
function. The above suggestions could be expanded greatly 
by changing the lattice in the assembly and repeating the 
particular investigation. Voids and absorbers could also 
be placed in the assembly to observe the effect of them upon 
the neutron wave . 
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Table 5. Original data at 0.63 radians/second (total 
number of counts) 



Oscillator 

position 

degrees 


One 

(3 min.) 


Two 

(3 sin.) 


Hole position 

Three tour 

(3 min.) (4 min.) 


Five 
(5 niirn ) 


0 


17,900 


8,623 


3,734 


2,254 


1,234 


30 


13,911 


8,711 


3,768 


2,220 


1,257 


60 


18,310 


8,600 


3,660 


2,170 


1,281 


90 


13,056 


8,224 


3,640 


2,105 


1,199 


120 


17*276 


8,189 


3,650 


2,102 


1,134 


150 


16,306 


7,84? 


3,522 


2,031 


1,074 


180 


16,026 


7,708 


3,374 


1,977 


1,042 


210 


15,995 


7,528 


3,368 


1,984 


1,100 


240 


16,248 


7,745 


3,484 


2,016 


1,077 


270 


17*005 


7,938 


3,491 


2,036 


1,036 


300 


16,971 


7,825 


3,496 


2,040 


1,103 


330 


17,747 


8,056 


3,605 


2,118 


1,164 



Table 6. 


Original data at 1.05 radians/second (total 
number of counts) 




Oscillator 




Hole position 




position 














One 


Two 


Three 


Four 


Five 


degrees 


(2 min.) 


(2 min.) 


(3 min.) 


(4 min. ) 


(5 min.) 


0 


11,950 


6,342 


3,963 


2,248 


1,297 


3° 


11,917 


6,036 


3,937 


2,291 


1,370 


60 


11,681 


5,721 


3,716 


2,165 


1,309 


90 


11,471 


5,366 


3,715 


2,215 


1,282 


120 


11,017 


5,351 


3,701 


2,193 


1,201 


150 


10,749 


5,230 


3,565 


2,019 


1,209 


180 


10,500 


5,223 


3,493 


1,968 


1,200 


210 


10,455 


5,198 


3,582 


1,984 


1,242 


240 


11,207 


5,233 


3,660 


2,034 


1,236 


270 


11,777 


5,28? 


3,554 


2,111 


1,268 


300 


11,218 


5,478 


3,769 


2,206 


1,229 


330 


11,373 


5,768 


3,879 


2,239 


1,249 
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Table 7. Original data at 1.63 radians/oecond (total 
number of counts ) 



Oscillator Hole position 

position 

One Two Three Four Five 

degrees (2 rain.) (2 rain.) (3 rain.) (4 rain.) (3 rain.) 



0 


12,179 


5,510 


30 


12,081 


5,513 


6o 


11,832 


5,426 


90 


11,471 


5,283 


120 


11,065 


5,282 


XfjG 


10,598 


5,114 


IGO 


10,625 


4,982 


210 


10,571 


5,000 


240 


10,946 


5,214 


270 


11,291 


5,191 


300 


11,665 


5,189 


330 


12,257 


5,234 



3,563 


2,082 


1,110 


3,479 


2,09? 


1,151 


3,448 


2,154 


1,096 


3,540 


2,067 


1,099 


3,444 


1,959 


1,113 


3,460 


1,932 


1,084 


3,310 


1,924 


1,076 


3,439 


1,894 


1,035 


3,369 


1,951 


1,028 


3,525 

3,479 


1,961 

1,356 

1,982 


1,039 

1,118 


3,499 


1,131 



Table 8. Original data at 2. 31 radians/second (total 
number of counts) 



Oscillator Hole position 

position 





One 


Two 


Three 


Four 


Five 


degrees 


(2 rain. ) 


(2 rain. ) 


(3 rain. ) 


(4 min. ) 


(5 rain. ) 


0 


12,027 


5,481 


3,656 


2,196 


1,235 


3° 


11,807 


5,407 


3,641 


2,110 


1,221 


60 


11,729 


5,428 


3,538 


2,177 


1,212 


90 


11,224 


5,354 


3,569 


2,048 


1,182 


120 


10,765 


5,361 


3,522 


2,034 


1,141 


150 


10,467 


5,197 


3,448 


2,000 


1,151 


180 


10,315 


5,153 


3,394 


2,002 


1,110 


210 


10,379 


5,252 


3,453 


1,973 


1,124 


240 


10,814 


5,298 


3,526 


2,092 


1,149 


270 


11,421 


5,333 


3,601 


2,150 


1,185 


300 


11,515 


5,424 


3,534 


2,136 


1,186 


330 


11,891 


5,436 


3,673 


2,194 


1,169 



72 



Table 9. Original data at 4.19 radians/ sec ond (total 
number of counts) 



Oscillator 

position 

degrees 


One 

(2 min. ) 


Hole position 

Two Three Four 

(2 min.) (2 min.) (4 min.) 


Five 
(5 min.) 


0 


11,590 


5,316 


2,417 


2,073 


1,124 


30 


IX, 585 


5,408 


2,396 


2,053 


1,145 


60 


11,376 


5,307 


2,319 


1,974 


1,142 


90 


10,095 


5,23b 


2,337 


1,967 


1,126 


120 


10,448 


5,326 


2,279 


1,989 


1,063 


160 


10,209 


4,975 


2,221 


1,944 


1,072 


180 


10,137 


4,857 


2,183 


1,884 


1,032 


210 


10,269 


4,896 


2,241 


1,949 


1,087 


240 


10,680 


4,827 


2,252 


1,952 


1,127 


270 


11,423 


4,992 


2,266 


2,001 


1,108 


300 


11,246 


5,222 


2,399 


1,994 


1,113 


330 


11,360 


5,406 


2,365 


2,103 


1,126 



Table 10. Original data at 6.28 radians/second (total 
number of counts ) 



Oscillator Hole position 

position 

One Two Three 



degrees (1 min.) (1 min.) (3 min.) 



0 


5,843 


3,211 


4,086 


30 


5,70S 


3,241 


4,051 


60 


5,680 


3,156 


4,012 


90 




3,120 


3,908 


120 


5,641 


3,013 


3,657 


150 


5,337 


2,989 


3,823 


ISO 


5,271 


2,981 


3,620 


210 


5,334 


3,050 s 


3,830 


240 


5,521 


3,129 


3,867 


270 


— 


3,130 


3,897 


300 


5,720 


3,188 


3,990 


330 


5,735 


3,245 


4,058 
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Table 11. Original data at hole position one (total number 
of counts) 



Oscillator 

position 

degrees 


, 7 ‘ 33 X 

(3 min.) 


Radians/sec ond 

9.42 12.55 

(2 min. ) (1 min. ) 


18.35 
{ 1 niin # ) 


0 


16,347 


11,261 


5,500 


5,367 


30 


16,474 


11,253 


5,452 


5,ol4 


60 


15*939 


10,843 


5,414 


5,359 


90 


15,665 


10,341 


5,354 


5,161 


120 


15,249 


10,458 


5,06? 


4,935 


150 


15,186 


9,962 


4,894 


4,879 


160 


15,208 


9,942 


b,9ol 


4,700 


210 


15,320 


10,140 


4,929 


4,916 


240 


15,932 


10,729 


5,148 


5,241 


270 


16,402 


10,970 


5,218 


5,326 


300 


16,390 


11,133 


5,508 


5,516 


330 


16,655 


11,196 


5,605 


5,308 



Table 12. 


Original data at hole position two (total number 
of counts) 


Oscillator 

position 

degrees 


, 7 - 33 ^ 
(3 sin. ) 


Radians/second 

9.42° 9-42 b 

(2 min. ) (2 min. ) 


12.55 
(1 rain.) 


0 


7,514 


5,18? 


5,182 


4,951 


30 


7,689 


5,058 


3,197 


4,906 


60 


7,602 


5,135 


5,0o| 


5,021 


90 


7,382 


5,042 


5,053 


%,9 78 


120 


7,298 


4,963 


5,110 


4,647 


150 


7,398 


4,975 


5,1H 


4,684 


180 


7,577 


4,829 


4,960 


4,573 


210 


7,368 


4 , 870 


4,972 


4,751 


240 


7,583 


4,975 


5,163 


4,750 


270 


7,677 


5,147 


5,239 


4,985 


300 


7,551 


5,332 


5,261 


4,93o 


330 


7,539 


5,232 


5,270 


5,063 



a First run. 
^Second run. 
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Table 13. Original data for calibration (counts/minute) 



Oscillator 














position 




Radians/sec ond 






degrees 


0*42 


0-63 


1.05 


1.68 


2.51 


. ‘i.ii 


0 


16,638 


19,215 


18,936 


18,700 


18,218 


18,003 


30 


17,374 


16,291 


18,338 


13,757 


18,269 


17,513 


6-0 


16,174 


17,270 


17,366 


17,475 


17,377 


17,016 


90 


14,729 


16,880 


16,140 


15,565 


15,615 


14,897 


120 


13,227 


13,609 


13,153 


12,384 


12,576 


12,442 


150 


9,155 


9,543 


9,701 


9,580 


9,453 


8,957 


160 


8,049 


7*947 


7,890 


7,790 


7,677 


7,460 


210 


8,573 


8,878 


8,651 


8,207 


8,371 


8,169 


240 


10,693 


11,497 


11,333 


10,989 


10,955 


10,726 


270 




14,412 


14,657 


14,104 


13,381 


13,591 


300 




17,078 


16,562 


16,604 


16,343 


16,223 


330 


«*» 


18,382 


18,198 


17,777 


17,763 


17,149 



degrees 


6*23 


7.33 


9.42 


12.55 


18.83 


0 


17,894 


17,897 


17,353 


17,542 


17,162 


30 


16,354 


17,781 


16,059 


17,128 


16,689 


60 


15,783 


16,361 


15,450 


15,812 


15,848 


90 


13,244 


14,208 


12,912 


13,748 


13,823 


120 


10,517 


11,001 


10,019 


10,790 


10,726 


150 


7,783 


8,505 


7,407 




7,923 


180 


7,171 


7,095 


6,939 


6,960 


6,709 


210 


8,419 


7,758 


8,401 


8,068 


8,144 


240 


11,594 


10,659 


11,391 


10,719 


10,844 


270 


14,907 


14, 066 


14,041 


13,696 


13,412 


300 


16,832 


16,303 


16,077 


15,682 


15,389 


330 


17,849 


17,642 


17,303 


16,967 


16,353 



